Dynamics of the bacterial flagellar motor with multiple stators 
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Abstract 



u 

The bacterial flagellar motor drives the rotation of flagellar filaments and enables many 
species of bacteria to swim. Torque is generated by interaction of stator units, anchored to 
the peptidoglycan cell wall, with the rotor. Recent experiments [Yuan, J. & Berg, H. C. 
(2008) PNAS 105, 1182-1185] show that near zero load the speed of the motor is independent 
^ ; of the number of stators. Here, we introduce a mathematical model of the motor dynamics 
that explains this behavior based on a general assumption that the stepping rate of a stator 



depends on the torque exerted by the stator on the rotor. We find that the motor dynamics 
can be characterized by two time scales: the moving-time interval for the mechanical rotation 
of the rotor and the waiting-time interval determined by the chemical transitions of the 
stators. We show that these two time scales depend differently on the load, and that their 



crossover provides the microscopic explanation for the existence of two regimes in the torque- 
speed curves observed experimentally. We also analyze the speed fluctuation for a single 
motor using our model. We show that the motion is smoothed by having more stator units. 
However, the mechanism for such fluctuation reduction is different depending on the load. 
We predict that the speed fluctuation is determined by the number of steps per revolution 
only at low load and is controlled by external noise for high load. Our model can be 
generalized to study other molecular motor systems with multiple power-generating units. 
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The swimming motion of bacterium Escherichia coli is propelled by the concerted ro- 
tational motion of its flagellar filaments [l|, 0]. Each filament (~ 10/im long) is driven by 
a rotatory motor embedded in the cell wall, with a angular speed of the order of 100 Hz 
0. The motor has one rotor and multiple stators in a circular ring-like structure roughly 
45nm in diameter [§| . The stators are attached to the rigid peptidoglycan cell wall and the 
spinning of the rotor drives the flagellar filament through a short hook (see [i] for a 3D 
reconstruction and Fig. 1(a) for a 2D sketch of the rotor-stator spatial arrangement). The 
rotor is composed of a ring of ~ 26 FliG proteins and each stator has four copies of proteins 
MotA and two copies of proteins MotB, forming two proton-conducting transmembrane 
channels. A flow of protons (or, in some alkalophilic and marine Vibrio species of bacteria, 
sodium ions), due to electrochemical gradients across the channels, causes conformational 
changes of the stator proteins that generate force on the rotor through electrostatic interac- 
tion between MotA and protein FliG The work per unit charge that a proton can do in 
crossing the cytoplasmic membrane through the proton channel is called the "proton-motive 
force" (pmf). 

At any given time, a stator is engaged with one of the 26 FliG monomers on the FliG ring 
as the duty ratio of the flagellar motor is close to unity 0. Presumably, the passage of protons 
switches the stator to be engaged with the next FliG monomer on the FliG ring along the 
direction of rotation, stretching the link between the stator and the rotor. The subsequent 
relaxation process rotates the rotor and the attached load towards the new equilibrium 
position. This can give rise to a step-like motion, characterized by advances of the rotor 
followed by waiting periods. The molecular details of the flagellar motor has been the 
subject of intense research^ and the step-like motion was recently demonstrated by direct 
observation 0] for a sodium-powered motor at very low pmf but a general understanding 
of the stepping dynamics of a single flagellar motor is still lacking. 

The torque-speed dependence is the key characteristics of the motor 0, Q, The mea- 
sured torque-speed curve (see Supporting Information (SI)) for bacterial flagellar motor 
shows two distinctive regimes. From its maximum value r max at stall (zero angular veloc- 
ity), the torque first falls slowly (by roughly 10%) as angular velocity increases at up to 
a large fraction (« 60%) of the maximum velocity, forming a plateau in the torque-speed 
curve. Then the torque starts to decrease quickly with increasing angular velocity, eventually 
approaches zero at the maximum velocity. For E. coli at room temperature under physio- 
logically relevant conditions, the maximum angular velocity is ~ 300Hz and the estimated 



maximum torque ranges from 4700 pN-nmllOjto 1400 pN-nm [ill, |12 



A few mathematical models [H, 14, 15.ll6l. 17. 18| have been proposed to explain various 
aspects of the observed torque-speed characteristics based on assumptions about details of 
the electrostatic interaction between the stators and the rotor. In a more general approach, 
recent work by Xing et al. [l9| has sought to understand the mechanism for the torque-speed 
curve characteristics without assuming a detailed description of the energy-transduction 
process. Their model can reproduce the observed torque-speed curve characteristics, and a 
set of general conditions to explain the observed torque-speed characteristics were suggested. 
However, the model by Xing et al. does not exhibit the correct behavior at low load. In their 
model, the maximum velocities depend inversely on the number of stators (see Supporting 
Text in [19]), whereas a recent experiment [2(| shows that near zero load the velocity of the 
motor is independent of the number of stators. 

Here, we aim at understanding both the torque-speed relationship and the individual 
motor dynamics by using a simple model describing the rotor's mechanical motion and the 
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stator's stepping probability. In our model, the stepping rate of a stator depends on the force 
between the stator and the rotor, in analogy to the Huxley model for Myosin (2l|. Specifically, 
"negative" force between a stator and its attached FliG monomer in the direction opposite 
to the rotation of the motor leads to a larger stepping rate for the stator. Under this 
general assumption, we find that the maximum velocity at (near) zero load in our model is 
mostly determined by the maximum stepping rate, independent of the number of stators, 
in agreement with the recent experiment by Yuan and Berg[2(|. Microscopically, the motor 
dynamics follows a repeated moving and waiting pattern characterized by two time scales: 
the moving-time interval t m associated with the (mechanical) rotation of the rotor and the 
waiting-time interval t w determined by the (chemical) transition of the stator. We find that 
t m and t w depend differently on the load and their crossover provides a natural explanation 
for the observed two regimes of the torque-speed curve. The fluctuation of the motor rotation 
are also studied in our model. We show that the sources of the motor speed fluctuation are 
totally different in the high and low load regimes and that the number of steps per revolution 
can only be extracted from the analysis of motor speed fluctuation in the low load limit. 



I. MODEL 

In Fig. 1(a), a schematic representation of a flagellar motor (rotor and stators) is shown. 
Each stator has two force-generating subunits symbolized by the light-blue and the red 
springs. The two force units of a stator interact with the FliG ring (rotor) in a hand-over- 
hand fashion as illustrated in Fig. 1(b), analogous to the way kinesin proteins interact with 



microtubules [22|, |23j. The switching of hands (force-generating unit) represents the energy- 
assisted transition when one hand releases its attachment and the other hand establishes 
its interaction with the FliG ring (rotor). The forces between the FliG ring and the stators 
drive the rotation of the rotor. In Fig. 1(c), the corresponding sequence of this hand- 
over-hand motion is shown in the energy landscape. The physical motion (solid arrow) of 
the rotor (green circle) is governed by its interaction potential with the engaged FliG. The 
hand-switch transition (dotted arrow) corresponds to a shift of the potential energy in the 
direction of motor rotation by angle So and the subsequent motor motion is governed by 
this new potential until the next switch. Microscopically, the shift angle could be different 
for the front hand and the back hand (with respect to the direction of the motor rotation); 
here for simplicity So is a constant. Commensurate with the periodicity S = 27r/26 of the 
FliG ring, we should have 2Sq = mS with a small integer m. In this paper we choose m — 1 
for simplicity. 

Due to the small Reynolds number, the dynamics of the rotor angle 9 and the load angle 
0l are over-damped and can be described by the following Langevin equations: 

^Tt = "lE V W ~ e *) ~ F (° ~ W + V^ff R a(t), (1) 



Ct-^ = F(0 - 9 L ) + V2k^ff L /3(t), (2) 

where £r and £l are the drag coefficients for the rotor and the load respectively, and N is 
the total number of stators in the motor. V is the interaction potential between the rotor 
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and the stator. V depends on the relative angular coordinates AOi = — Of, where Of is 
the internal coordinate of the stator i. Of increases by Sq when the stator switches hands. 
This discrete change in Of is called a jump of the stator in this paper. The load is coupled 
to the rotor via a nonlinear spring described by a function F, which can be determined 
from the hook spring compliance measurement of ref. [24!] (see Fig. SI in SI). The last 
terms in Eq. ([TJ12]) are stochastic forces acting on the rotor and on the load, with fc# the 
Boltzmann constant, T the absolute temperature, and a(t) and f3(t) independent white noise 
fluctuations of unity intensity. 

The dynamics of the stator i is governed by the transition probability for the discrete 
jump of its internal variable Of during the time interval t to t + At: Pi(0f — > Of + 5q). In this 
paper, P, is assumed to depend on the torque generated by the z'th stator n = — V'(AOi), 
which depends on the relative angle AOf 

Pi{0 S i -> Of + So) = r(r t )At = k(A0i)At. (3) 

The specific form of the jumping rate r(rj) (or k(AOi)) is unknown. We assume it to be 
a decreasing function of 7$, with the stator stepping rate being higher when is negative 
("Tj < 0, AOi > 0) than when is positive (7$ > 0, AOi < 0). 

Fig. 1(d) illustrates the motor dynamics, where the rotor (green circle) is either pulled 
forward or dragged backward by individual stators (purple circles) depending on their rela- 
tive coordinates with respect to the rotor. The stator coordinate changes by jumping forward 
by So with a probability rate that is a function of its relative coordinate. For simplicity, we 
set the potential function V to be a V^-shaped function: V(A0) = to\A0\, and the torque 
from a single stator is tq with its sign depending on whether the stator is pulling (AO < 0) 
or dragging (AO > 0). Correspondingly, the stator jumping rate depends on the sign of the 
force: k(A0 < -S c ) = 0, k(-S c < AO < 0) = k + , k(A0 > 0) = £;_(> k +) as illustrated in 
Fig. 1(e). A cutoff angle S c is introduced to prevent run- away stators. Quantitatively, we 
use r = 505pN-nm, £ R = 0.02pN-nm-s-rad" 1 , k + = 12000s" 1 , k_ = 2k + , S c = S in this 
paper unless otherwise stated. The load £l varies from 0.002 — 50pN-nm-s-rad _1 . Simulation 
time step At = 0.01 — lfis. 

II. RESULTS 

A. Two characteristic time scales and their different dependence on the motor 
speed 

In Fig. 2(a), a typical case of time dependence of the rotor angle 0(t) from our model 
is shown. The motion of the rotor consists of two alternating phases: moving and waiting. 
The moving phase occurs when the net force on the motor is positive (in the direction of 
motion). The waiting phase is when the system reaches mechanical equilibrium (net force 
equals zero) and the motions are driven by thermal fluctuation. The dynamics of the motor 
can thus be characterized by the two time scales t m and t w . The waiting-time t w is the 
time the rotor spends fluctuating around a equilibrium position, i.e., the bottom of the total 
potential Vt = V(A0i). Once in the waiting phase, the rotor can only start to have a 
net motion when a stator jumps to break the force balance and thus shift the equilibrium 
position forward. The subsequent net motion of the rotor to reach the new equilibrium 
position takes t m , which is defined as the moving-time. The definitions of t w and t m are 
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shown in Fig. 2(b). 

The dynamics of the motor depend on the load, higher load leading to slower speed. We 
study how the two scales t m and t w vary with the load or equivalently the speed of the 
motor (speed is chosen because of its direct measurability in experiments). We find that 
the two time intervals have very different dependence on the motor speed as shown in Fig. 
2(c). The waiting-time interval is determined by independent chemical transitions, i.e, by a 
Poisson process with rate k, so we have: (t w ) oc (k^ 1 ). Since k varies between two constants 
k + and A;_ (except for extreme high load where k — 0), the averaged waiting-time has only 
a weak dependence on motor speed as shown in Fig. 2(c). On the other hand, the average 
moving-time can be estimated as: (t m ) « 5 m /uj m , with 5 m the average angular movement, 
uj m = r m /(fa + fa) the average speed, and r m = (— V/) m the average net torque in the 
moving phase. Increasing the load fa leads to a decrease of the speed uj m and an increase 
of the moving-time. In addition, at lower speed, it is more likely for stator to jump in the 
middle of a moving phase before the system reaches its force equilibrium. These premature 
stator jumps effectively increase 5 m and further increase (t m ). These two factors lead to a 
strong dependence of (t m ) on the load (or the speed) as shown in Fig. 2(c). Besides the 
difference in their average values, the distribution functions for t m and t w are also different 
(See Fig. S2 in SI for details). 



B. The two regimes of the torque-speed curve 

In Fig. 3, the torque-speed curves calculated from our model for 8 different stator numbers 
are shown. Our model results closely resemble the observed torque-speed curves. There is 
a plateau regime with almost constant (10% decrease) torque from zero up to a large speed 
(^ 100Hz for N = 8), followed by a steep declining regime of the torque, all the way to zero 
at a speed of roughly 300Hz. By using the two time scales t m , t w , and noting that the net 
torque is zero during the waiting phase of the motor, the time-averaged torque r and speed 
uj can be estimated: 

(tm) $m i , \ 

T ~ ~, T ~, r T m, U « ; ; r. (4) 

(t m ) + (O (tm) + (O 

The two distinctive regimes in the torque-speed curve can be understood intuitively within 
our model by the different dependence of (t m ) and (t w ) on the speed shown in the last 
section. 

In the low-speed (high-load) regime defined by (t m ) ^> (t w ), we have r r m and uj ~ 
b m /(t m ) from Eq.(|3]). As discussed in the last section, for low speed a stator can jump 
prematurely during the moving phase before the system reaches the bottom of the potential 
well. As a result, each stator spends most of its time generating positive torque r . Therefore, 
in this high-load regime, while the speed changes significantly, the torque stays near its 
maximum value T max = Nto, which is proportional to the number of stators. 

In the high-speed (low-load) regime defined by (t m ) (t w ), we have r w r m (t m ) / (t w ) 
and uj ~ 5 m / (t w ) from Eq.(jl]). As shown in Fig. 2(c), for increasing speed (t m ) decreases 
quickly while (t w ) remains roughly the same. This naturally explains the steep decrease of 
the torque r with speed in the high-speed regime. Intuitively, in this high-speed regime, a 
stator can be pushed into the negative torque region (A9 > 0) because the rotor rotates 
too fast for the premature jump to occur. As the stators spend large fractions of their time 
dragging the rotor, the torque of the motor decreases quickly. 
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The different dependence of the waiting and moving-time intervals on the speed not 
only gives a clear general explanation for the two regimes of the torque-speed curve, it also 
explains the sharpness of the transition between the two regimes. Since the dependence of 
(t m ) on the speed is much steeper than that of (t w ) (as shown in Fig. 2(c)), the crossover 
between the two regimes takes place in a small region of the speed values, thus making 
the two regimes in the torque-speed curve well defined, as found in both experiments and 
simulations of our model. 



C. Independence of the motor speed on the number of stators at near zero load 

At near zero load, our model shows that the motor moves with a roughly constant speed 
that is independent of the number of stators, as demonstrated in Fig. 3. Recent res- 
urrection experiments using gold nano-particle (extremely low load) indeed showed such 
independence [20)]. The mechanism for this surprising behavior can be understood with our 
model. In the low-load regime, the motor spends most its time in the waiting phase where 
the net torque is zero. In our model with symmetric potential V, this force equilibrium 
is achieved by having on average half of the stators pulling the rotor and the other half 
dragging it. If we number the stators in Fig. 1(d) from left to right, the rotor's equilibrium 
position sits between the iV/2'th and the (N/2 + l)'th stators. This equilibrium state breaks 
under two possible scenarios: 1) One of the N/2 dragging stators jumps to the pulling side. 
This occurs with a probability rate Nk_/2. 2) The (N/2 + l)'th stator jumps and shifts 
the equilibrium to a position between the iV/2'th and the (N/2 + 2)'th stator. This occurs 
with probability rate k + . The average distance between the new and the old equilibrium 
positions are <5 m (~ Sq/N) and 5 m /2 for scenario 1) and 2) respectively. The fundamental 
reason for the decrease in step size with N is due to the high duty ratio as first recognized 
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25j, |26|. Similar step size reduction with N was recently observed in kinesin-1 motor [27 



The maximum speed oo max near zero load is then: 

fcW(JV) « S m ^ + 5 fk + » [1 + (k + /k_)N-% (5) 

which only depends weakly on N, if k + /k_ <C 1. The estimated maximum speed 0J max oc 
k_5 /2 makes sense as u) max should be limited by the step size and the maximum stepping 
frequency of an individual stator. 

We have studied the dependence of u max on the ratio r = k + /k^ and iV by numerical 
simulations of our model. In Fig. 4(a), we show the torque-speed curves for N — 1 and 
N = 8 for two different values of r: r = 0.2 and r = 1.2. To quantify the dependence 
of u max on N, we define a quantity A = 2(u max (l) - uj max (8)) / (u max (l) + u max (8)) to 
characterize the relative difference between the maximum speeds for motors with one and 
eight stators. As shown in Fig. 4(b), u max is roughly independent of N, i.e., |A| < 0.1 as 
long as r < 0.5. However, u: max (l) becomes significantly bigger than uj max (8) for r > 1. The 
observed dependence of A on r agrees well with the analytical estimate given by Eq.©. 
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D. Motor speed fluctuation at different load levels and the estimate of step num- 
bers 



The measured motor speed fluctuates due to two main factors: the external noise such as 
the Brownian noise and measurement noise, and the intrinsic probabilistic stepping dynamics 
of the stators. Samuel and Berg[26| first investigated the speed fluctuations by studying the 
smoothness of the periodic motor motion characterized by T = n(Ti) 2 /((T 2 ) — (T n ) 2 ), where 
T n is the period for n revolutions. By measuring T in a resurrection experiment where 
the stator number is inferred from the discrete increments in average motor speed, it was 
found that T is proportional to the number of stators. The proportionality constant was 
interpreted as the number of steps per revolution. Here, we analyze the motor fluctuation 
by using our model to understand how different noise sources contribute to T and how T 
behaves differently at different load levels. 

For low loads, the motor spend most of its time in the waiting phase. The average motor 
step size is 5 m ~ <5 /A^(<^ 27r), there are n s = 2tx /5 m ~ 2txN/8q steps in each revolution, and 
the average periodicity is (Ti) = n s (t w ). Since the waiting-time intervals are uncorrected, 
the variance of the n— revolution periodicity can be expressed as: (T 2 ) — (T n ) 2 = nn s ((t^) — 
(t w ) 2 ). Furthermore, because the waiting-time t w is determined by a Poisson process, its 
variance is equal to (t w ) 2 . F can thus be written as: 



(T 2 ) - (T n ) 2 (tl) - (t w ) 2 5 

showing that T = ^yN is proportional to the stator number N, and the proportionality 
constant 7 = 2h/5q corresponds to the number of steps per stator per revolution. This 
behavior is verified in our model by calculating T during a simulated resurrection process, 
where additional stators are added by a Poisson process with time constant t s = 400s (Fig. 
5). For near zero load, the average speed is independent of the stator number in agreement 
with [20[ (see Fig. 5(a)). However, T increases by a fixed amount 7 = 2tt/5 as a new stator 
is incorporated into the system as shown in Fig. 5(b), consistent with the analytical result 
by Eq. ([()]). The behavior of T as shown in Fig. 5(b) represents a quantitative prediction of 
our model that could be tested in resurrection experiments with extreme low load, such as 
in 0- 

For high load, the net torque is roughly constant r ps Ntq and the speed can be expressed 
as ojq = t/(£l + Cr) ~ Ntq/^l, which explains the constant increment of speed for every 
additional stator (up to eight) seen in our model (Fig. 5(c)) as well as in the resurrection 
experiments by Blair and Berg[28|. For additional stators beyond a certain large number 
of stators, the speed passes the knee in the speed-torque curve and our model predicts a 
decrease in the speed increment, which is consistent with the recent experiments by Reid 



et al [12| that showed the same decrease in speed increment as stator number goes up 
to AT = 11. The dynamics of the load angle can be obtained by summing Eqs.(l-2) and 
taking the limit £r/£l — » 0. This leads to: 9l — u>o + ^/2K B T/£ L [3(t), which describes 
the simple motion of the load with a constant speed t^o perturbed by random noise. From 
the equation for 9l, the periodicity and its variance can be determined: (T n ) w 2rni/ij , 
(T 2 ) — (T n ) 2 ps 4mi k bT / (^lUq). We can now express T as: 

r = ^ K^COgK^N (7) 

L -(T 2 )-(T n ) 2 kgT k B T iy - [<) 
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r is again proportional to N through its dependence on the speed loq. However, unlike in 
the low-load regime, the proportionality constant 7 = tttq/ (ksT) has nothing to do with the 
number of steps per revolution. Instead, 7 depends on the ratio between the intrinsic driving 
force (to) and the external noise ksT, as the driving force overcomes the external noise to 
make the motor moves smoothly. This behavior is verified in our model by calculating T 
during the resurrection simulation. As shown in Fig. 5(d), T goes up with the number of 
stators but with a much larger proportionality constant 7, which quantitatively agrees with 
the expression 7rr /(kBT) from our analysis. 

Therefore, although the motor-speed fluctuation is always suppressed by higher numbers 
of stators, the mechanisms are different for different load levels. For low load, the smoother 
motion for larger N is caused by the increase in step number per revolution. For high load, 
the smoother motion for larger iV is caused by larger driving force (therefore larger speed) 
in comparison with the constant external noise. The difference in motor fluctuation between 
the high and the low load regimes is confirmed by our simulation as shown in Fig. 6, where 
the proportional constant T is shown for different values of external noise strength /c#T (Fig. 
6(a)) and different load (Fig. 6(b)). 



III. SUMMARY AND DISCUSSION 



We have presented a mathematical description of the rotary flagellar motor driven by 
hand-over-hand power-thrusts of multiple stators attached to the motor. All key observed 
flagellar motor properties [2], including those from a recent resurrection experiment at near 
zero load[20j, can be explained consistently within our model. The crucial ingredient of 
our model is that the hand- switching rate depends on the force between rotor and stator. 



This feature is known to be valid for other molecular motors, including kinesin[29| and 
myosin [30j] . Therefore our model should be generally applicable to the study of these linear 
motors, especially in the case when there are multiple power-generating units attached to 
the same track [2 71]. 

For the flagellar motor, we find that its dynamics follows an alternating moving and 
waiting pattern characterized by two time scales t m and t w . The mechanism underlying 
the observed torque-speed relationship and its dependence on the number of stators, is 
revealed by studying the dependence of these two time scales on the load. For high load, 
(t w ) <S (tm), the motor spend most time moving (albeit slowly) with all the stators pulling 
the motor in the same direction. So the torques generated by individual stators are additive, 
leading to a roughly constant torque r max ~ Nr , which persists up to the knee speed u n . 
Microscopically, the existence of this torque-plateau regime is due to the premature stator 
jumps which prevent the stators from going into the negative torque region. Since the rate of 
the premature jumps is k(A9 < 0), larger k + and larger cutoff S c increase the knee-speed u n 
(see Fig. S3 in SI for details). For low load, (t w ) ^> (t m ), the motor spends most time in the 
waiting state. A waiting period ends when one of the dragging stators jumps to the pulling 
side or the pulling stator closest to the bottom of the potential well jumps. Therefore the 
maximum motor speed io max is limited by the maximum jumping rate of the stators. uj max 
can be estimated from our model. Eq.© shows that u max has only a weak dependence on 
N for small k + /k_, as confirmed by simulations of our model, and in agreement with recent 
resurrection experiments at near zero load [20]. Eq.([5]) also explains the strong dependence 



rT 

of 

oJmax on A in a recent model by Xing et al[19|. The jumping probability used in [1 
has a complicated profile and is maximum in the positive torque region. In our model, this 
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would correspond to having k + /k^ ^> 1, which is the opposite to what is required to achieve 
independence of u max on N. 

The robustness of our results were verified using different forms of the rotor-stator po- 
tential V and the force function F between the load and the rotor. In particular, we have 
studied a smoothed symmetric potential with a parabolic bottom and an asymmetric po- 
tential V (similar to the one used in [l9|) where the negative torque (r_) is bigger than the 
positive torque (t + ). We find that all of our general results remain the same (see SI and Fig. 
S4&S5 for details). For the asymmetric potential, the condition for u max being independent 
of N is generalized to k+/k- <C t + /t_ . From the analysis and direct simulation of our model, 
we do not find any significant dependence of the torque-speed curve characteristics on the 
specifics of the force function F between the load and the rotor. In particular, contrary to 



what was proposed in [19J], there is no difference between the case of a viscous load that 
interact with the rotor through a soft spring and a viscous load without spring (see Fig. S6 
in SI for details). 

Besides the torque-speed curve which describes the time averaged behavior of the motor, 
we have also studied the speed fluctuation for individual flagellar motor. We find that the 
fluctuation is damped by the number of stators for all load levels. However, we show that 
the dominating source of the motor fluctuation is different depending on the load. For low 
load, the speed fluctuation is dominated by the discrete stochastic stepping events whereas 
for the high load, it is controlled by the external noise, such as Brownian fluctuations or 
possibly measurement noises. The original measurements on motor fluctuation by Samuel 
and Berg [2^] were done in the high- load regime as evidenced by the discrete increment of 
speed in their resurrection experiment. Therefore the strength of the fluctuations obtained 
there is probably more reflective of the strength of external noise than the number of steps 
per revolution. It would be interesting to perform the fluctuation analysis in the low-load 
regime as achieved in 2(| to determine the steps number per revolution and compare with 
the recent direct observation of the steps [(J. 

Simple relations between the macroscopic observables ( 

Tmaxi ^ni ^max) &nd the micro- 
scopic variables of the system (to, k+, fc_) are established by analysis of our model. These 
relations can be used to predict the microscopic parameters quantitatively from the torque- 
speed measurements. They can also be used to study the dependence of the flagellar motor 
properties on other relevant external parameters such as the pmf the temperature, and 
solvent isotope effects. For example, since changing of pmf gives rise to self-similar torque- 
speed curves|32], we conclude from our model that larger pmf not only increases the chem- 
ical transition rates fe's, it also increases the stator- rotor interaction strength tq. Changing 
temperature or replacing H + with D + (solvent isotope effect) should affect the chemical 
transition rates. These changes in k + and fc_ lead to changes in the knee speed u n and the 
maximum speed uj max in our model without changingthe maximum torque at stall, which 



is consistent with previous experimental observations [9l. 1311 . 

Backward stator jumps with 6 s — > 9 s — 5 can be incorporated in our model to study 
the relatively rare motor back-steps [g]. The back-jumps are neglected in this paper as their 
probabilities are much smaller than those for the forward jump in the region of relative 
angles (Ad > —5 C ) relevant for our study here. However, we expect the back-jumps to 
become dominant for A8 < — 5 C where the forward jumps are prohibited. Since the landing 
points of these back-jumps are still on the positive side of the potential with positive torque 
tq, inclusion of back-jumps in our model for A9 < —5 C can naturally explain the observed 



torque continuity near stall when the motor is driven backwards by an optical tweezer 10 
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In our model, the step size depends inversely on the stator number N. This behavior 
is a general consequence of duty ratio being unity and independent stepping of the stators, 
as pointed out by Samuel and Bergjl^]. This A -1 dependence of the step size seems to be 
inconsistent with an "apparent independence" of step size on N claimed in [6|]. However, a 
careful study of the experimental data reveals that the N dependence of the step size can 
not be ruled out, because the step size distribution was measured for a varying population of 
stators, whose number was neither controlled nor measured precisely in [6|. An unambiguous 
way to determine whether the step size depends on N is to measure the step size for different 
fixed N or at least to measure N simultaneously. Such experiment was done recently for 



kinesin-l[27|] and showed that step size for iV = 2 is half of that for N = 1. 

Our model works for the clockwise (CW) as well as for the counterclockwise (CCW) 
rotation. It was recently suggested that the switching between the CW and CCW state of 
the motor is a non-equilibrium process and the energy needed to drive the motor switch 



could be provided by the same pmf that drives the mechanical motion of the motor [33 
The possible link between the switching pro cess and the rotational motion of the motor is 
supported by experimental observations [11] showing that the average switching frequency 
depends on the proton flux. It is therefore highly desirable to develop an integrated model 
to describe both the mechanical part of the flagellar motor, associated with the rotational 
motion, with the signaling part, associated with the switching process. More experimental 
information on the components of the motor (M-ring/C-ring/MotAB) and how they interact 
with each other are needed to achieve this goal. 

IV. ACKNOWLEDGMENTS 

This work is partially supported by a NSF grant (CCF-0635134) to YT. 



10 



[1] Berg, H. C. and Anderson, R. A. (1973) Bacteria swim by rotating their flagellar filaments, 
Nature 245, 380-382. 

[2] Berg, H. C. (2003) The rotatory motor of bacterial flagella. Annu. Rev. Biochem. 72, 19-54. 
[3] Thomas, D. R. Francis, N. R. Xu, C. and DeRosier, D. J. (2006) The Three-Dimensional 

Structure of the Flagellar Rotor from a Clockwise-Locked Mutant of Salmonella enterica 

Serovar Typhimurium. J. Bacteriol. 188, 7039-7048. 
[4] Kojima, S. and Blair, D. F. (2006) Conformational change in the stator of of the bacterial 

motor. Biochemestry 40, 13041-13050. 
[5] Ryu, W. S. Berry, R. M. and Berg, H. C. (2000) Torque-generating units of the flagellar motor 

of Escherichia coli have a high duty ratio. Nature 403, 444-447. 
[6] Sowa, Y. Rowe, A. D. Leake, M. C. Yakushi, T. Homma, M. Ishijima, A. and Berry, R. 

M. (2005) Direct observation of steps in rotation of the bacterial flagellar motor. Nature 437, 

916-919. 

[7] Sowa, Y. Hotta, H. Homma, M. and Ishijima, A. (2003) Torque-speed Relationship of the 

Na+-driven Flagellar Motor of Vibrio alginolyticus. J. Mol. Biol. 327, 1043-1051. 
[8] Chen, X. C. and Berg, H. C. (2000) Torque-Speed Relationship of the Flagellar Rotatory 

Motor of Escherichia coli. Biophysical J. 78, 1036-1041. 
[9] Berg, H. C. and Turner, L. (1993) Torque generated by the flagellar motor of Escherichia coli. 

Biophysical J. 65, 2201-2216. 
[10] Berry, R. M. and Berg, H. C. (1997) Absence of a barrier to backwards rotation of the 

bacterial flagellar motor demonstrated with optical tweezers. Proc. Nat. Acad. Sci. USA 94, 

14433-14437. 

[11] Fahrner, K. Ryu, W. S. and Berg, H. C. (2003) Bacterial flagellar switching under load. Nature 
423, 938. 

[12] Reid, S. W. Leake, M. C. Chandler, J. H. Lo, C-J. Armitage, J. P. and Berry, R. (2006) The 
maximum number of torque-generating units in the flagellar motor of Escherichia coli is at 
least 11. Proc. Nat. Acad. Sci. USA 103, 8066-8071. 

[13] Lauger, L. (1988) Torque and rotation rate of the bacterial flagellar motor. Biophys. J. 53, 
53-65. 

[14] Meister, M. Caplan, S. R. and Berg, H. C. (1989) Dynamics of a tightly coupled mechanism 
for flagellar rotation. Bacterial motility, chemiosmotic coupling, protonmotive force. Biophys. 
J. 55, 905-914. 

[15] Berry, R. M. (1993) Torque and switching in the bacterial flagellar motor. An electrostatic 

model. Biophys. J. 64, 961-973. 
[16] Braun, T. F. Poulson, S. Gully, J. B. Empey, C. Van Way, S. Putnam, A. and Blair, D. 

(1999) Function of Proline Residues of Mot A in Torque Generation by the Flagellar Motor of 

Escherichia coli. J. Bacteriol. 181, 3542-3551. 
[17] Walz, D. and Caplan, S. R. (2000) An electrostatic mechanism closely reproducing observed 

behavior in the bacterial flagellar motor. Biophys. J. 78, 626-651. 
[18] Schmitt, R. (2003) Helix Rotation Model of the Flagellar Rotary Motor. Biophys. J. 85, 

843-852. 

[19] Xing, J. Bai, F. Berry, R. and Oster, G. (2006) Torque-speed relationship of the bacterial 
flagellar motor. Proc. Nat. Acad. Sci. USA 103, 1260-1265. 



11 



Yuan, J. and Berg, H. C. (2008) Resurrection of the flagellar rotatory motor near zero load, 
Proc. Nat. Acad. Sci. USA 105, 1182-1185. 

Huxley, A. F. (1957) Muscle structure and theories of contraction. Prog. Biophys. Biophys. 
Chem. 7, 255318. 

Asbury, C. L. Fehr, A. N. and Block, S. M. (2003) Kinesin Moves by Asymmetric Hand-Over- 
Hand Mechanism. Science 302, 2130-2134. 

Yildiz, A. Tomishige, M. Vale, R. D. and Selvin, P. R. (2004) Kinesin walks hand-over-hand. 
Science 303, 676-678. 

Block, S. M. Blair, D. and Berg, H. C. (1989) Compliance of bacterial flagella measured with 
optical tweezers. Nature 338, 514-517. 

Samuel, A. D. T. and Berg, H. C. (1995) Fluctuation analysis of rotational speeds of the 
bacterial flagellar motor. Proc. Nat. Acad. Sci. USA 92, 3502-3506. 

Samuel, A. D. T. and Berg, H. C. (1996) Torque-Generating Units of the Bacterial Flagellar 
Motor Step Independently. Biophysical J. 71, 918-923. 

Leduc, C. Ruhnow, F., Howard, J. and Diaz, S. (2007) Detection of fractional steps in cargo 
movement by the collective operation of kincsin-1 motors. Proc. Nat. Acad. Sci. USA 104, 
10847-10852. 

Blair, D. F. and Berg, H. C. (1988) Restoration of torque in defective flagellar motors. Science 
242, 1678-1681. 

Visscher, K.Schnitzer, M. J. Block, S. M. (1999) Single kinesin molecules studied with a 
molecular force clamp. Nature 400, 184-189. 

Oguchi, Y. Mikhailenko, S. V. Ohki, T. Olivares, A. O. De La Cruz, E. M. and Ishiwata, S. 
(2007) Load-dependent ADP binding to myosins V and VI: Implications for subunit coordi- 
nation and function. Proc. Nat. Acad. Sci. USA 105, 7714-7719. 

Chen, X. C. and Berg, H. C. (2000) Solvent-Isotope and pH Effects on Flagellar Rotation in 
Escherichia coli. Biophysical J. 78, 2280-2284. 

Gabel, C. and Berg, H. C. (2003) The speed of the flagellar rotary motor of Escherichia coli 
varies linearly with protonmotive force. Proc. Nat. Acad. Sci. USA 100, 8748-8751. 
Tu, Y. (2008) The nonequilibrium mechanism for ultrasensitivity in a biological switch: Sens- 
ing by the Maxwell's demons. Proc. Nat. Acad. Sci. USA 105, 11737-11741. 



12 



StatQT 




o 1 2 eis 




Stators that drag Stators that pull 



FIG. 1: A Model for the flagellar rotary motor, (a) Schematic illustration of the rotor-stators 
spatial arrangement. The rotor contains 26 FliG proteins and there are multiple stators, each with 
two subunits (red and light blue springs), (b) A sequence of three rotor-stator configurations (from 
top to bottom) illustrating the hand-over-hand interaction between the two subunits of a stator 
and the FliG proteins in the rotor, (c) The same sequence as in (b) is shown in the potential 
landscape. The solid arrow represents the physical rotation of the rotor angle (9) down a given 
(V-shaped) potential, the dotted arrow represents the chemical change (switching of hands) that 
shifts the potential, (d) The full motor model with multiple stators in the angle space. Each stator 
is represented by its internal angle 6 s . The rotor is pulled forward by the stators in front it and 
dragged back by the stators behind it. The stator angle can only change by jumping forward with 
rate k that depends on the relative angle A9 = — s . The form of k(A6) used in this paper is 
given in (e), which shows the dragging stators have a higher jump rate k- > &;+ and a cutoff angle 
-5 C where k(A6 < -5 C ) = 0. 
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(a) 




FIG. 2: The motor dynamics and its dependence on the load, (a) Rotor angle 9 versus time 
for ~ 4/5 of a revolution for N = 1. The angular unit is the FliG periodicity 5. The insert is 
enlarged in (b). (b) Zoom of the time series in (a) showing two complete steps. Solid line shows 
the stator position 6 s . A jump in 6 s marks the start of a moving phase for the rotor and the 
waiting phase starts when the rotor catches up with the stator. The definitions of the moving-time 
t m and the waiting-time t w are shown, (c) The average waiting-time (t w ) (dashed line) and the 
average moving-time (t m ) (solid line) over 500 revolution as a function of the rotational speed for 
N= 1. (t w ) decreases slowly with increasing speed from 1ms to 0.1ms while (t m ) decreases much 
faster from roughly 50ms to 0.005ms. 
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FIG. 3: The torque-speed (r — u) curves for different stator numbers (N = 1 to N = 8) from our 
model. Two regimes of the r — to curves, i.e., constant r up to a large knee speed oj n and fast 
decrease of r to zero at the maximum speed to ma x, ar e evident for all stator numbers. The torque 
per stator t/N versus the speed oj is shown in the insert. The torque at stall scales with N while 
the maximum speed is independent of N. 
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FIG. 4: Dependence of the maximum speed at zero load on N as a function of r = k+/k-. (a) 
Torque-Speed curves for r = 0.2 (solid lines) and r = 1.2 (dotted lines) are shown for N = 1 (black 
lines) and N = 8 (red lines). We change r by varying k + and keeping k- constant. For r = 0.2, 
the maximum velocities ui max (1) and 

^maxi^) 3>t zero load are roughly the same, while they differ 
significantly for r = 1.2. (b) Xhe dependence of A = 2(w ma:r (l) ^max (8))/(uw(l) + (8)) 
on r. The red line represents our analytical predictions from Eq.(5). The shaded region shows the 
~ 12% experimental error. 
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FIG. 5: The motor speed and its fluctuation in a simulation of the resurrection process for low 
and high loads, (a) Speeds as a function of time after successive stators are added at low load 
(£l = 0.002pN-nm-s-rad -1 ). The simulation shows no dependence of the speed on the number of 
stator (labeled with a number form 1 to 8). (b) The smoothness parameter T with n= 5 calculated 
from the time series shown in (a). T value is calculated with a moving-time window of 1.5s (~ 500 
revolutions), (c) Same as (a) at high load = 8pN-nm-s-rad _1 ) showing the roughly linear 
dependence of the motor speed u> on N. (d) V at high load from the time series shown in (c) using 
a window of roughly 10s corresponding to a total of 100 revolutions, (b) and (d) show that V 
increases with N at both high and low load. 
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FIG. 6: The speed fluctuation and its dependence on the load (with N = 1). (a) The dependence 
of T on the external noise strength defined as T/Tq, where To is the room temperature. T depends 
strongly on T for high load (crosses) while it is a constant determined by the step size at load load 
(dots), in agreement with our theoretical results (lines). (b) The dependence of T on speed oj for 
T = Tq. Typical error bar (SD) is shown at high load; the error at low load is comparable to the 
symbol size. 
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